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The LEP Collaborations have reported a small excess 
of events in their combined Higgs boson analysis at center 
of mass energies \fs < 208 GeV. In this communication, I 
present the result of a calculation of the probability distribu- 
tion function of the Higgs boson mass which can be rigorously 
obtained if the validity of the Standard Model is assumed. It 
arises from the combination of the most recent set of precision 
electroweak data and the current results of the Higgs searches 
at LEP 2. 

PACS numbers: 14.80. Bn, 12. 15. Mm, 12.15.Ji. 



Combining all Higgs decay channels and experiments, 
the LEP Collaborations report a 2.5a excess in their 
data H. The probability that this is due to an upward 
fluctuation of the background is 0.6%. Of course, the 
Higgs boson has been searched for in many different en- 
ergy bins, and there is an infinitely large energy range out 
of reach, so that one expects to observe an upward fluc- 
tuation somewhere. It is therefore difficult to interpret 
these numbers, and it would be imprudent to conclude 
that the Higgs boson has been found with 99.4% proba- 
bility. 

In this communication, I present the answer to a differ- 
ent but related question, which is Given the data, what 
is the probability that the Higgs boson is within reach of 
LEP 2? This question can be answered unambiguously 
once the probability distribution of the Higgs boson mass, 
Mh, has been constructed. This is not possible given the 
Higgs search results at LEP 2 by themselves, regardless of 
how strong a signal is observed there. The reason is that 
there is an infinite domain of a priori possible values of 
Mh beyond the kinematic reach of LEP 2. As a result, 
the Mh distribution is improper, i.e., it is asymptoti- 
cally non-zero. Including the electroweak precision data, 
however, renders it sufficiently convergent and a proper 
integration is possible. 

The electroweak precision data by itself, especially Z 
pole asymmetries and the W boson mass, is now precise 
enough to constrain Mh significantly, as can be seen from 
Fig. [I]. Including the latest updates as presented at the 
2000 summer conferences, I find from a global fit to all 
data using the package GAPP B , 
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Note, that by definition the central value in Eq. (|l|) maxi- 
mizes the likelihood, Ne~ x ( M «)/ 2 ; and that correlations 
with other parameters, , are accounted for, since mini- 
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FIG. 1. Constraints on Mh from various sets of precision 
data as a function of the top quark mass. For the individ- 
ual data sets I show lcr contours, while the ellipse for all 
data refers to the 90% CL. The excluded region by LEP 2, 
M H < 113.2 GeV (95% CL), is also indicated. 



mization w.r.t. these is understood, y 2 = Xmin- The error 
is the standard lcr uncertainty (A% 2 = 1). The 68.27% 
central confidence interval, 



53 GeV <M H < 131 GeV, 



(2) 



differs slightly from Eq. (|l|) due to the non-Gaussian 
(asymmetric) distribution of In My. The 90% central 
confidence interval yields, 



38 GeV < M H < 173 GeV. 



(3) 



The right hand side of Eq. (f|, i.e., the 95% upper limit, 
does not take into account the direct search results of 
the LEP Collaborations. Negative search results will in- 
crease the upper limit, because the probability distribu- 
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tion function is effectively being renormalized. For ex- 
ample, in a previous analysis || we found that the Higgs 
exclusion curve presented by the LEP Collaborations in- 
creased the 95% upper limit by 30 GeV. 

The use of the Higgs exclusion curve, however, is only 
appropriate if no indication of an excess is observed. In 
general, it is more appropriate to consider the likelihood 
ratio for the data, 

£(data|signal + background) 

Qlep 2 = rl , 7 n — ] t; 1 (4) 

£(data| background) 

where both the numerator and the denominator are func- 
tions of Mff. The quantity, 



AX 2 (direct) = —2 In Qlep 2 



(5) 



can then be added to the x 2 -function obtained from the 
precision data. If the signal hypothesis gives a better 
(worse) description of the data than the background only 
hypothesis we find a negative (positive) contribution to 
the total x 2 - Note, that this is a consistent treatment 
also in the case of a large downward fluctuation of the 
background or even if no events are observed at all. Use 
of Qlep 2 had been originally advocated in Ref. § (see 
Eq. (23) in that reference). 

This treatment can be rigorously justified within the 
framework of Bayesian statistics [^,0, which is particu- 
larly suited for parameter estimation. Bayesian methods 
are based on Bayes theorem M , 



p(Mff |data) 



p(data\M H )p(M H ) 
p(data) 
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which must be satisfied once the likelihood, p(data|Mjj), 
and prior distribution, p{Mn), are specified. p{data) = 
J p(data|Mff )p(Mff)<iMff in the denominator provides 
the proper normalization of the posterior distribution on 
the left hand side. Depending on the case at hand, the 
prior can 

1. contain additional information not included in the 
likelihood model, 

2. contain likelihood functions obtained from previous 
measurements, 

3. or be chosen non-informative. 

Of course, the posterior does not depend on how infor- 
mation is separated into the likelihood and the prior. As 
for the present case, I choose the informative prior, 



P(M H ) = Qlep 2 P non ~ inf (M H ), 



(7) 



where the non-informative part of the prior will be chosen 
as 



This choice corresponds to a flat prior in the variable 
In Mh, and there are various ways to justify it Q. One 
rationale is that a flat distribution is most natural for a 
variable defined over all the real numbers. This is the 
case for \\iMh but not M H . Also, it seems that a priori 
it is equally likely that Mh lies, say, between 30 and 
40 GeV, or between 300 and 400 GeV. In any case, the 
sensitivity of the posterior to the (non-informative) prior 
diminishes rapidly with the inclusion of more data. As 
discussed before, p(Mh) is an improper prior but the 
likelihood constructed from the precision measurements 
will provide a proper posterior. 

Occasionally, the Bayesian method is criticized for the 
need of a prior, which would introduce unnecessary sub- 
jectivity into the analysis. Indeed, care and good judge- 
ment is needed, but the same is true for the likelihood 
model, which has to be specified in any statistical model. 
Moreover, it is appreciated among Bayesian practition- 
ers, that the explicit presence of the prior can be advan- 
tageous: it manifests model assumptions and allows for 
sensitivity checks. From the theorem (^|) it is also clear 
that any other method must correspond, mathematically, 
to specific choices for the prior. Thus, Bayesian methods 
are more general and differ rather in attitude: by their 
strong emphasis on the entire posterior distribution and 
by their first principles setup. 

Including Qlep 2 in this way, one obtains the 95% 
CL upper limit Mh < 201 GeV, i.e. notwithstanding 
the observed excess events, the information provided by 
the Higgs searches at LEP 2 increase the upper limit by 
28 GeV. 

Given extra parameters, £ l , the distribution func- 
tion of Mh is defined as the marginal distribution, 
p(Mjj|data) = Jp(M H , f |data) U t P(€) d C- If the like- 
lihood factorizes, p(Mn, £*) = p(Mh)p(C)> the C depen- 
dence can be ignored. If not, but p(C\Mh) is (approxi- 
mately) multivariate normal, then 



X 2 (M ff ,f) = X 2 nin (M ff )+ 



i d 2 x 2 (M H ) 

2 d&dtj 



(C - & in (MHM j - din(^))- 



The latter applies to our case, where £ l 
(m tl a sl a(Mz)). Integration yields, 



p(Afe|data) ~ VdeT~Ee- x — 



(9) 



where the error matrix, E — ( ° . g)f g ^ ) \ introduces 
a correction factor with a mild Mh dependence. It cor- 
responds to a shift relative to the standard likelihood 
model, x 2 (M H ) = xlin( M H) + A X 2 (M H ), where 



A X 2 (M H ) = In 



det E(M H ) 



p non ~ ini (M H ) = M H \ 



(8) 



det£(M z ) ' 

This effect tightens the Mh upper limit by 1 GeV. 
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FIG. 2. Probability distribution function for the Higgs 
boson mass. The probability is shown for bin sizes of 1 GeV. 
Included are all available direct and indirect data. The shaded 
and unshaded regions each mark 50% probability. 



I also include theory uncertainties from uncalculated 
higher orders. This increases the upper limit by 5 GeV, 

M H < 205 GeV (95% CL). (11) 

The entire probability distribution is shown in Fig. |^. 
Taking the data at face value, there is (as expected) 
a significant peak around M H = 115 GeV, but more 
than half of the probability is for Higgs boson masses 
above the kinematic reach of LEP 2 (the median is at 
Mh = 119 GeV). However, if one would double the in- 
tegrated luminosity and assume that the results would 
be similar to the present ones, one would find most of 
the probability concentrated around the peak. A similar 
statement will apply to Run II of the Tevatron at a time 
when about 3 to 5 fb _1 of data have been collected. 

The described method is robust within the SM, but it 
should be cautioned that Mh extracted from the preci- 
sion data is strongly correlated with certain new physics 
parameters. Likewise, the Higgs searches at LEP 2 de- 
pend on the predictions of signal and background expec- 
tations which are strictly calculable only within a spec- 
ified theory. This note focussed on the Standard Model 
Higgs boson. 
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